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Abstract 

The photon production arising due to time variation of the medium has 
been considered. The Hamilton formalism for photons in time-variable medium 
(plasma) has been developed with application to inclusive photon production. 
The results have been used for calculation of the photon production in the course 
of transition from quark-gluon phase to hadronic phase in relativistic heavy ion 
collisions. The relative strength of the effect as well as specific two-photon cor- 
relations have been evaluated. It has been demonstrated that the opposite side 
two-photon correlations are indicative of the sharp transition from the quark- 
gluon phase to hadrons. 



1 Introduction 



The formation of the quark-gluon plasma (QGP) in relativistic heavy ion coUi- 
sions with subsequent transition to hadrons is under discussion for many years. 
Numerous measurable signals from QGP phase (such as J/V^ suppression) have 
been suggested. At the same time, the observation of enhanced antibaryon pro- 
duction together with lattice calculations which predict low temperature 
of the transition are indicative of the fast QGP-hadron transition 0] without for- 
mation of the mixed phase. In this paper we consider a new specific mechanism 
of photon production which is effective in the case of sharp transition from QGP 
to hadrons, see Q]. 

The phenomenon under consideration is the photon production in the course 
of evolution of strongly interacting matter. Let us consider photons in the medium 
at initial moment to having momentum k and energy Um. Let the properties of 
the medium (its dielectric penetrability) change within time interval St so that 
the final photon energy is Uf. As a result of the energy change the production of 
extra photons with momenta ±k takes place these photons having specific two- 
photon correlations. Moreover the photons with opposite momenta and opposite 
helicities are produced even in the absence of the initial photons. Analogous 
processes were considered for mesons [§], §|, 05 HI and applied to pion production 
in high-energy heavy ion collisions 0. The conditions for strong effect are the 
following: first, the ratio of the energies ujin/uf must not be too close to unity 
and second, the transition should be fast enough. 

The calculation of the effect requires consideration of the Hamilton equations 
of motion. So in Section 2 we present in short the corresponding formalism which 
is extended in Section 3 to the case of time-dependent medium. A simple method 
of calculation of the photon correlations is presented in Section 4 with subsequent 
estimation in Section 5 of the polarization operator and the evolution parameter 
which determines the the strength of the transition effect. In the last Section 6 
we calculate the transition effects in heavy ion collisions and present the results 
of the work. 

2 Basic formulation 

We are interested in time evolution of photon creation and annihilation operators 
a|(k, t) and aj(k, t) in the medium. In this section the medium is taken in its rest 
frame. The properties of the medium (plasma) will be described by the transverse 
polarization operator n(ct;, k, T, /i, m) which depends on energy u;, momentum k, 
temperature T = 1/(3, chemical potential and the mass m of the charged 
particles of the plasma. To get the evolution law of a|(k, t) and aj(k, t) we will 
use the Hamilton formalism. 

The Lagrangian of the electromagnetic field in the medium is taken in the 
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form 

L^\jd\ (e(x, i)eE(x, t) - tf(x, t)) (1) 

where dielectric penetrability e acts as the factor e(a;, k) in the momentum space 
(otherwise it acts as an operator). The object of quantization is the real valued 
vector potential A(x, t). In our case we can use the gauge conditions 

ylo(x,t) = 0, d«t;A(x,t) = (2) 

so that 

E(x, t) = -A(x, t), H(x, t) = rotA(x, t) (3) 

The second of Eqs.(2) means that the vector potential is transverse one (kA = 0) 
having two components Ai (usual linear polarizations). The transverse dielec- 
tric penetrability e is connected with photon transverse polarization operator 11 
through equation 

e{uj,k) ^l-U{uj,k)/u'^ (4) 

where n(a;, k) will be cosidered as the real valued and even function of cu and k 
(see Section 5). Turning to momentum representation 

A(x,t) = J^^e-^'°'+^^-A{k,ko) (5) 



and using (1-5) the action of the system can be written in the form 

S^lJ d^kM-k, -ko) [kl - k^ - n(k, ko)] A{k, ko) (6) 

Variation of the action provides the equation 

'kl - k^ - n(k, ko)] A{k, ko) = 0, 1 = 1,2 (7) 

and the energy cu of the photon in the medium is given by the usual dispersion 
equation 

-k^ -U{LU,k) = (8) 

Let us note in advance that Il{uj, k) will be positive and slowly varying function 
of k playing essentially the role of the effective photon mass squared. 

Introducing the time and momentum dependent field coordinates q(k, t) 

and coming one step back in (6) we get the Lagrangian in the form 
L = ^Jd^k [qi{-k, t)qi{k, t) - k\i{-k, t)qi{k, t) - qi{-k, t)fl{k)qi{k, t)] (10) 
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where the polarization operator acts either to the left or to the right giving 
equivalent results. In this form the polarization term is related to potential 
energy (unlike q terms) ensuring the physically sensible form of the Hamiltonian 
in the case under consideration. The Lagrange equations are given by variation 
of (10). Evidently they have the oscillator form 

g,(k)+u;2^,(k) = (11) 

with u determined from (8). That is 

fl(k)g,(k) = n(^,k)g,(k) (12) 

for plane wave solutions of the equations of motion. 

Turning to Hamilton formalism for quantum fields, we introduce the canoni- 
cally conjugated momentum 

Oqi{K t) 

and postulate canonical equal time commutation relations 

fe(ki,i),p,(k2,i)] = iM(ki - ks) (14) 

with all other commutators being zero. Let us note that the presence of the 
opposite sign of k in the rhs of (13) is necessary also for compatibility of (14) 
with commutation relations in coordinate space 

A(xi,t),ij(x2,t) 

Below this opposite sign of the momentum k will result in important physical 
consequences describing production of photon pairs with opposite directions of 
momenta of the photons. 

The Hamiltonian is introduced in the usual way 

j d^kpi{k)qi{k) - L 
^IJd'k (p,(-k)p,(k) + k2g,(-k)g,(k) + g,(-k)fl(k)g,(k)) (16) 

and the Hamilton equations for Heisenberg operators 

q,{k)=t[H,q,ik)]=p,{-k), 
Pi{-k) = I [H,pi{-k)] = -uj^iik) (17) 

connect qiik) withpi(— k) being consistent with (11), (12). 



i5ij5{yii - X2) 



(15) 
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Let us at last introduce the photon creation and annihilation operators a'^{k.) 
and a(k) through decomposition 

(2^(2^ («.(k,t)e^^- + aKk,t)e~^^-) (18) 

They are defined here for stationary medium in initial and final states with con- 
stant u) and they are connected in these regions with canonical coordinates and 
momenta by equations 



:(-k)=zy|(al(-k)-a,(k)) (19) 



as it can be seen from comparison of representations (9) and (18) and their time 
derivatives. The creation and annihilation operators have simple time dependence 

a{t) ~ e~"^\ a\t) ~ e^"^* (20) 
corresponding to running waves (photons) and satisfy commutation relations 

a,(ki),a](k2)] = 5,,(5(ki - ka) (21) 
which follow from canonical commutation relations (14). 

3 Photon evolution in time dependent 
medium 

Let the polarization operator be the time dependent function of its parameters. 
In this case the Hamilton equations (17) remain valid having time dependent 
energy a; (t) (as well as Lagrange equation (11) which is their consequence). This is 
confirmed by the fact that their solution represents the time dependent canonical 
transformation conserving commutator (14). Indeed, the solution of (17) can be 
written in the form (cf | 



Qiik, t) = siit)qiik, 0) + S2it)pii-k, 0) 
Qiik, t) = pi(-k, t) = si(t)g,(k, 0) + S2{t)pi{-k, 0) (22) 

where Si{t) and S2{t) are two linearly independent real valued solutions of the 
classical equation (11) with time dependent energy uj and with initial conditions 

si(0) = l, ii(0) = 0, S2(0) = 0, S2(0) = l (23) 
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One can see from (22) that the canonical commutator transforms in the following 
way: 

[g,(k,t),p,(k,t)] = Wisi,S2) [g.(k,0),p,(k,0)] (24) 

where 

W{s,, S2) = s,{t)s2it) - S2(t)ii(t) (25) 

is the Wronskian determinant of (11) which does not depend on time for this 
equation (due to absence of the term with the first order derivative in the equa- 
tion), this constant determinant being equal 1 due to initial conditions(23). 

Let us consider the process of evolution of the photons from the initial state 
with energy uji to asymptotic final state with energy uj2. The final state annihi- 
lation and creation operators are given by 



ai(k,t) 



{u2qi{k,t) +ipi{-k, t)) 



al(k, t) = --^= {LJ2qi{-k, t) - ipi(k, t)) 



(26) 



as it follows from (19). We substitute solutions (22) for qi(k,t),pi{k,t) and 
introduce two lineary independent complex valued classical solutions C,(t),^*(t) 
instead of Si{t), S2(t) 

e(t) = si{t) + iuJiS2{t), C{t) = Slit) - liuMt) (27) 

with initial conditions 

e(0) = r(0) = l, ^(0) = z^i, C{t) = -iooi (28) 

Then, using (19) for the initial state, we get the Bogoliubov transformation |]Tl| 
connecting the creation and annihilation operators in the initial and final states 
(let us remind that these operators were defined only for asymptotic states having 
constant energy u): 



ttiik, t) = u{k, t)ai{k, 0) + v(k, t)aj(-k, 0) 
al(k, t) = v*ik, t)ai{-k, 0) + u*{k, t)al{k, 0) 



(29) 



with 



em + -{•(() 

UJ2 



1 

2' 



v{k, t) 

It follows from (30) that 

u*{t)u{t) -v*{t)v{t) 



m + -m 



2uj- 



(30) 



(31) 
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where 

w{^,c)^mt{t)-im*{t) (32) 

is again the time independent Wronskian determinant. So 

u*u -v*v = l (33) 

due to initial conditions (28). In turn, as it follows from (29), the condition (33) 
ensures conservation of the commutator 

aiik, t), at(k, t)] = [a,(k, 0), a|(k, 0)] (34) 

(the last commutator is equal to V/{27r)^ for the system having volume V). In 
view of condition (33) the Bogoliubov coefficients u, v can be represented in the 
form 

u = cosh r(k)e^"i , v = sinh r (k)e^"2 (35) 

where r(k) is the main parameter which determines the photon production and 
the phases ai,a2 do not play important role and they will not be considered 
below. 

To find the coefficients u, v (for fixed momentum k) one must turn to classical 
equations for oscillator having variable frequency (energy). We look for solution 
of the equation 

^ + u\t)^^0 (36) 
In the initial state, where the energy cui is constant, we take a single wave 

m = e'^'\ t<Un = (37) 

At large enough time, t > tf, when the energy 002 becoms constant again, the 
general solution has the form 

^(t) = + Cae-^'"^*, t > tf (38) 

Substituting (37), (38) into (30) we get corresponding Bogoliubov coefficients 

u = flcie-'''^\ V = J—C2e-''^^\ t > tf (39) 

V ^1 V '^i 

where constants Ci, C2 satisfy relationship 

I c, \' - 1 C2 r= ^ (40) 

UJ2 

due to condition (33). 

To find the final expressions for coefficients u, v one must know the full solution 
of (36) connecting the asymptotics (37) and (38). Here one can use an analogy 
between the above problem and the problem of the wave propagation through 
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(above) the one-dimensional potential barrier. In the last case (38) (after sub- 
stitution of X for t) represents incoming and reflected waves and (37) represents 
outgoing wave. The reflection coefficient 



C2/C1 1 2= I v/u 



tanh r 



(41) 



gives the desired ratio of the Bogoliubov coefficients (up to phases). Therefore 
one can use known quantum-mechanical results. So, shifting the initial time tin 
to large negative value, we take the reference model of the energy variation: 



(42) 



The problem with such form of the potential barrier can be found in text- 
books |T^. It contains the important parameter 8t giving characteristic time 
of the energy variation. The evolution parameter r is given now by 



tanh ^ \'vl'u\ 



1 /tanh(7ru;25r/4)\ 

2 ^ I tanh(7ru;i5r/4) / 



(43) 



For sharp transition [iOibr <^ 1) we get from (43): 



In 



UJ2 



UJ2 — UJi 
UJ2 + UJi 



UiSr <^ 1 



(44) 



This is the case of the most intensive pair production. For large transition time 
6t the process becoms adiabatic one and the evolution parameter falls down 
exponentially: 



-T7UJ2St/2 



-ituii5t/2 



uj6t > 1 



(45) 



Let us note that the above results (44) for sharp transition can be found imme- 
diately from Hamilton equations (17). Suggesting that the canonical coordinate 
q and momentum p are finite, we see from these equations that q and p are finite 
as well and so q and p are continuous functions of time at the transition point: 



q{-6t) = q{6t), p{-5t) = p{6t), 6t^0 



(46) 



coinsiding for both sides of the sharp border between media with different pho- 
ton energies uJi and uj2. Being expressed according (19) through creation and 
annihilation operators a^, a at both sides of the border, the equations (46) lead to 
Bogoliubov transformation (29) between pairs of operators ,a at different sides 
of the border with coefficients 



\u\ 



UJ2 + OOi 



2y/UJiUJ2 

The last equations correspond to (44). 



UJ2 - UJl 



2y/UJiUJ2 



(47) 
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4 Inclusive photon production in a simple 
model 



Below we will be interested in inclusive photon production in heahy ion colli- 
sions. We confine ourselves to symmetric case when the photons having opposite 
momenta ±k are produced in an equivalent way (central collisions of identical 
nuclei). For simplicity the Bogoliubov coefficients u{'k),v(k.) will be taken to be 
real valued and /c = |k| dependent. To get feeling of the main features of the pho- 
ton production and their correlations (and for further references and comparison) 
we formulate in this section a simple model - fast simultaneous transition of large 
homogeneous system at rest (the movement of the system will be considered in 
Section 6). 

For physical interpretation of the evolution effect it is helpful to introduce the 
complex valued vectors of the photon circular polarization: 

e± = (eiiies)/ A ke± = (48) 

and corresponding components of operators a(k), a^(k) 

a± = (oi ±ia2)/V2, ai ^ {a{ ^ ial) / V2 (49) 

so that 

a(k) = e+a_(k) + e_a+(k), a± = e±a (50) 

The components a±, a± satisfy the standard commutation relations of the form 
(21) with i = ± and represent the creation and annihilation operators of the pho- 
tons with definite spin projection ±1 on the direction of the photon momentum 
(the helicity). 

In what follows we will denote the creation and annihilation operators in the 
final state as b"^ ,b leaving notations a^, a for the operators in the initial state. In 
terms of the components b± the helicity operator is 

Sbs = iib\b2 - bib,) = blb+ - blb^ = Nt,+ - N,_ (51) 

for every momentum k, the photon number operator is 

Nb = blb+ + bib- = Nb+ + Nb- (52) 

and the Bogoliubov transformation takes the form: 

b±{k) = ua±{k) + vai{-k), 4(k) = u*ai{k) + v*a±{-k) (53) 

where indexes ± are helicities related to corresponding momenta k or — k. 

The resulting average number of the photons and their correlations depends 
on the initial state. We suggest that the initial state is the statistical (gaussian) 
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mixture of the coherent states |T3|, |5[ of the photons of each polarization. Then 
using (53) we get the photon momentum distributions 

= (iV,i(k)) = (l + \v\') (iV,±(k)) + \v\'{K^{-k)) + \v\'-^^ (54) 



As one can see from (54) the photons are created in pairs having opposite 
directions of their momenta and opposite spin directions (the same hehcities). If 
both kinds of polarizations are equally represented in the initial state for every 
momentum k, then the same property is valid for the final state, 

{SUk)) = (55) 

and the photons produced with opposite momenta have equal helicities. The in- 
tensity of the transition production is given by the Bogoliubov coefficient |f p and 
the last term in rhs of (54) represents the result of the ground state rearrangement 
( the initial ground state is not the ground state for 6-operators which operate 
in another medium in the final state). We suggest that k — k symmetry takes 
place. Then the resulting photon momentum distribution can be written in the 
form: 

^ = (6i(k)6,(k)) = ^ [n{k) + {2n{k) + 1) sinh^ r] (56) 

(sum over helicities A) where ?T,(k) is the average level occupation number of the 
single mode in the initial state and we used the parametrization (35). The photon 
production in the course of the transition is given by the second term in rhs of 
(56). It is weak for small evolution parameter r{k) being of the order r^. 

The evolution effect is better seen in photon correlations. Two-photon inclu- 
sive cross-section is given by 

^ (6l(ki)6Uk2)&A(ki)6,(k,)) = (6l(k06A(ki))(6t (k2)6,(k2)) 



+ (6l(ki)6,(k2))(6t(k2)6A(ki)) + (6i(k06i(k2))(6A(ki)6,(k2)) (57) 

The first term in rhs of (57) is the product of single-photon distributions, the sec- 
ond term gives the Hanbury Brown-Twiss effect (HBT, called also Bose-Einstein 
correlations) and the third term is essential if the time evolution effect takes place 
giving opposite side photon correlations. 

Production of the photons with opposite directions of their momenta was 
already noted after (54). However it was implicitely suggested there that the 
volume of the system is arbitrary large, V ^ oo. In real situation we deal with 
large but finite size of the colliding nuclei. The finite size of the photon source 
will smooth out the effect. The same is valid for HBT effect responsible for 
the same side correlations. The special technique was elaborated in a number 
of papers to describe the source size effect for HBT correlations including the 
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use of Wigner phase space density [0] and the method of equivalent classical 
currents [|15| which can be used here. Instead we directly modify the creation 
and annihilation operators in such a way that they are nonzero only inside 
some region in coordinate space which is described by the function /(x): 



6a(x) = 6a(x)/(x), 6i(x) = 6l(x)/(x) (58) 
Then the modified (smoothed out) operators in momentum space are: 

6i(ki) = I d'kfik^ - k)6l(k), h{k,) = I d'kf{k - k2)6A(k) (59) 

where the Fourier transform /(ki — k) is sharply peaked around the point k = kj 
(smoothed 5-function). 

The introduced operators satisfy modified commutation relations 



6A(ki),6Uk2 

d^x 
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Xv 



d3fc/(ki-k)/(k-k2) 



5xu 



J-^jW''^-'''^^- = 6,Mk^-k, 



(60) 



Here F{ki — k2) is the form-factor of the photon source (Fourier transform of 
the source density) and it also represents the smoothed 5-function which has the 
width of the order of the inverse size of the source. In particular 



F(0) = V/{2n) 



(61) 



where V is the effective volume of the source at the stage of the transition in 
the medium. Let us note that above we in fact already used (61) substituting it 
for 5^(0) in the contribution of the ground state rearrangement to the evolution 
effect in (54). This factor appears also when one wants to use the level population 
function n(k) (see (56)), having its origin in the correspondence of discrete and 
continous Fourier decompositions. Using operators modified in three-dimensional 
coordinate space we suggest that the volume V changes inessentially in the course 
of the time transition. 

Now one can apply the Bogoliubov transformation (53) to estimate the cor- 
relators of the modified operators: 



(6^(ki)&±(k2)) = / d^k 



u\k)n{k) + v\k){n{k) + 1)] /(ki - k)/(k - k2) 

^ n(k) + (2n(k) + 1) sinh^ rl ^(ki - k2) (62) 



and in a similar way 



V 



(6±(ki)6±(k2)) = (6^(ki)&i(k2)) 
■{2n{k) + 1) sinhr(A;) coshr(A;)G'(ki + k2) 



(63) 
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with other correlators vanishing in the course of statistical averaging in the case of 
sharply peaked form-factor. In the above equations G(ki ± is the normalized 
form-factor (G(0) = 1) and we took into account that the functions f(ki — k) 
and /(k — k2) as well as form-factor G(ki ± k2) are sharply peaked functions 
of their arguments (at zero momentum) having characteristic scale of the order 
of inverse size of the source, this scale being much less than the characteristic 
scales of the momentum distribution n{k) and the evolution parameter r{k). So 
the last two functions can be evaluated at any of momenta ki,k2 ~ ±k (we 
suggest that the process is k — >• k symmetric). The equations (62-63) together 
with (57) show that the correlations arising due to photon identity (HBT effect) 
are the same side momentum correlations of the photons having the same helic- 
ities whereas the photons arising due to transition effect have the opposite side 
momentum correlations and approximately opposite spin directions (the same 
helicities again). 

Returning to two-photon correlations given by (57) (sum over helicities) we 
get the relative correlation function which is measured in experiment: 

C(ki, k2) = 1 + ^G\k^ - k2) + ^Rl{k)G\k, + k2) (64) 

with 

(2n(A;) + 1) sinhr(A;) coshr(A;) 
°^ ' ~ n{k) + {2n{k) + l)sinhV(A;) ^ ^ 

according to (62-63). As can be seen from (64-65) the transition effect depends 
strongly on the evolution parameter r(k),see(43-45). Below, after necessary mod- 
ifications, we apply the above considerations to photon transition radiation in 
heavy ion collisions. 



5 Photons in plasma and the evolution 
parameter 

To find the evolution parameter r{k) one must know the photon energy in plasma.The 
spectrum of photons in plasma is given by dispersion equation 

col = k'^ + Il{uJk,k,T, fx,m) (66) 

Here 11 is the polarization operator for transverse photons dependent on temper- 
ature T = chemical potential n and the mass m of charged particles. Below 
we use an approximate form of 11 extracted from original expression : 

U{uj,k)=ujl 



1 - 



k^ 



UJ 



vk 



(67) 
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with 

,2_ 4(7« r^^f^2 ^2o2V/^ 



dx (^x^ — w? jS'^^ nF{x,fiP) (68) 



where a = 1/137, v"^ is the averaged velocity squared of the charged particles in 
the plasma, factor g takes into account the number of the particle kinds and their 
electric charges {g = 5/3 for u,d quarks) and np is the occupation number of the 
charged particles: 

npiM = (e^--^^ + l) ' + (e^"+^'^ + l)"' (69) 

(Fermi distribution). The polarization operator for scalar charged particles is ap- 
proximately a half of that for fermions with substitution of Bose distribution for 
Fermi distribution. The imaginary part of the polarization is small in the case 
under consideration (no Landau damping) and it was neglected in the calcula- 
tions. 

For small masses of charged particles,/3m <^ 1, the asymptotics of the polar- 
ization operator takes the known simple form 

This form can be used in quark-gluon phase. For large masses of fermions (con- 
stituent quarks and nucleons) when /?m > l(m > /x) it is convenient to use the 
expansion 

^l(F) = ^ E cosh(n/3/i)ifi(n/3m) (71) 

n=l ^ 



where Ki{x) is the exponentially decreasing modified Bessel function and the 
main contribution comes from a few first terms of the series. Analogous expansion 
was also used for bosons (pions) having no chemical potential: 



2 



Aam ^ 1 



J2-K,inl3m) (72) 



We considered the polarization operator and photon spectrum for three pos- 
sible kinds of plasma: quark-gluon plasma (QGP) with u, d light quarks, con- 
stituent quark(m = 350MeV^)-pion plasma and hadronic (pions and nucleons) 
plasma. Chemical potential (baryonic one) was taken to be equal to lOOMeV 
per quark corresponding to typical value for CERN-SPS energies (say IGOGeV per 
nucleon in Pb — Pb collisions). The temperature Tc of the transition was taken 



to be equal to rather small value lAOMeV characteristic for final hadrons |]T7 
(Tc = 200MeV was also tested giving close results). Under above conditions the 
asymptotic values Ua in (70-72) are equal to 24, 11.7, 11.2 and 2.5 in MeV units 
for light quarks, constituent quarks, pions and nucleons correspondingly. 
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Let us note that the approximation (67) where we use the averaged velocity 
squared v"^ = (v^), suggests that higher moments of the velocity distribution do 
not differ significantly from corresponding powers of f^. So we calculated the 
ratio {v'^) / {v'^Y verified that it differs from unity inessentially being 1 for 
QGP, 1.11 for constituent quarks (valons), 1.06 for pions and reaching 1.25 for 
nucleons which last contribution is small by itself. Corresponding values of v'^ 
used in our estimations are 1, 0.545, 0.731 and 0.300. 

Evidently the polarization operator in (66) plays the role of (momentum de- 
pendent) photon mass squared m^. The termal mass squared at zero photon 
momentum k is equal to 

= ^l{^- y) (73) 

and it approaches uj^ at large momenta. The slope of the dispersion curve at the 
origin is: 



varying from 1 at w = to 1.2 at v = 1. In rough approximation the polarization 
operator taken at the dispersion curve can be represented by the simple expression 

n(u;,,A;) = m^(^)=u;o' + ^|^, ^ = ^ " ^ ' ^^^^ 

reproducing position of the point cjq) asymtotic value uja and the slope (74). In the 
presence of two different sorts of charged particles (say pions and nucleons) the 
expression (75) should be modified in an appropriate way to take into account the 
presence of two contributions and to ensure the new correct slope doJ^ jdk^. The 
simple approximation (75) appears to be rather close to polarization operator (67) 
taken at the dispersion curve (66) (these two curves have a common point also at 
some finite A;) and it will be used below for estimation of the evolution parameter 
r{k). Zero momentum photon masses m^(0) was found to be 19.5, 14.4 and 10.0 
in MeV units for QGP, valon-pion and nucleon-pion plasma correspondingly. 

The evolution parameter r(A;), which gives the strength of the transition radi- 
ation, was determined in (43-45) for the reference model (42). It depends strongly 
on transition time 5r. We suggest that this time interval is not large being of the 
order of 1/m/c. Then for small momenta k {u!{k)dT <^ 1) the parameter r{k) is 
universal and it can be well approximated by the expression which follows from 
(44): 

_ m^,(fc)-m^,(fc) _ ^m^(fc) 
^^^^ " 4((m2(A;)) + k') ~ 4{u;'{k)) ' ""^"'^^ ^ ^ 

where m^j are photon termal masses squared at both sides of the transition 
and (^^) is their average. Zero momentum evolution parameter r(0) is equal to 
0.324, 0.154 and 0.178 for QGP-hadron, QGP-valon and valon-hadron transitions 
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correspondingly. Higher momentum behaviour of r{k) is model dependent. Below 
it will be taken in the form 

r(A;) = ^exp(-|Mr) (77) 

The Eq.(77) is a simple version of the asymptotical form (45) and it can be 
well sewed together with (76) giving a monotonically decreasing function of the 
momentum. Let us note that on general grounds one expects that r{k) falls 
down exponentially at large kSr if the time dependence of the energy uj{k,t) in 
the course of the transition has no singularities (non-analiticity) at real time axis. 
Below Eqs. (76-77) will be used for estimation of the transition effect in heavy ion 
collisions. Only QGP-hadron transition will be considered. In view of fact that 
the evolution parameter r{k) appears to be small number at all momenta fc, all 
expressios will be estimated in the lowest order in r{k). 



6 Transition effect in heavy ion collisions 

Let us now apply the above consideratins to photon production in heavy ion 
collisions. Let us suggest that the quark-gluon plasma is formed at the initial 
stage of the ion collision. Let the plasma undergoes expansion and cooling. The 
expansion is taken to be longitudinal and boost invariant []TB|. Recent lattice 
calculations show rather low critical temperature of the deconfinement and 
chiral phase transition, ~ ISOMeV^ as well as sharp drop of the pressure up to 
very small value when the temperature approachs T^. thus provocing instability 
in the presence of overcooling. So we do not expect long-living mixed phase and 
consider fast transition from quark to hadron matter with characteristic transition 
proper time duration 5t of the order of Ifm/c. 

To calculate the transition effect one must shift to rest frame of each moving 
element of the system and integrate over proper times r and space-time rapidities 
1] of the elements of the system. Then the invariant single-photon distribution in 
central rapidity region y = reads: 



dN 



(PkTdy 



J^QGP + -'fa 



Hr 



rdr / dr] 



d'MPo^) +JdvJ d'xT^,r\p) (78) 



d^p 



with p = kr cosh t]. 

The first term in rhs of (78) describes photon production from hot quark- 
gluon plasma. Here is the QGP production rate per unit four-volume in the 



rest frame of the matter 19 
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= 18^ exp(-p/T) ln(l + — ) (79) 

with OL = 1/137, as = 0.4, k = 0,58. It can be used also for hadron gas as 
its uncertainity is larger than the difference between the first-order QGP and 



hadron gas production rates . Contribution from hadronic resonances are not 
considered here. The second term in rhs of (78) describes photon production due 
to transition from QGP to hadrons in the vicinity of proper time Tc,cf (56). The 
time duration of the transition is taken to be small in this term in comparison 
with total time duration of photon production process. The evolution parameter 
r(p) is given here by (76-77). 

As the last step one must specify the temperature evolution in (78). We 
suggest that the temperature depends on proper time of the volume element 
with power-like dependence: 

(T/To) = (r/ro)-'/' (80) 

where tq and Tq are initial proper time and initial temperature. For final estima- 
tion we use 6 = 3 typical for hydrodynamical picture and choose low transition 
temperature Tc = 140MeV. After transition the photons live some time in 
hadronic medium and we suggest termal momentum distribution of the hadrons 
(modified by the expansion of the system). We neglect termal photon production 
below Tc and do not introduce a special freese-out temperature. An alterna- 
tive way is to introduce the final temperature Tf < Tc. Not pretending to high 
accuracy we shall not distinguish between and Tf below. 

Using now the (79-80) the photon production rate in QGP phase in (78) can 
be integrated over space-time rapidities rj: 

W^|;;^.,2.)V^.J|^/.V.«c. (81) 

dx- -TTTT \n{l + Kx-\ (82) 

The remaining integral (82) over temperatures can be easily evaluated in different 
subregions of the photon momentum kx- Below we will be interested in rather 
low photon momenta kx (up to 500Mef) where transition effect is expected 
to be more pronounced. In this momentum region the integral (82) depends 
mainly on final temperature Tc (at asymptotically large momenta kx it depends 
on initial temperature Tq). The final proper time in (81) depends on initial 
conditions in general. However for two main variants of the initial conditions used 
in literature l2l[|:roTo = l,To/Tc = 3/2(ro = 1/m/c) and tqTo = 1/3, Tq/Tc = 
5/2(ro = 0.2/m/c) the proper time Tc changes inessentially, being 3.125/m/c and 



with 

'QGP 
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3.375/m/c correspondingly (for 6 = 3). Below we will use for an average value 
Tc = 3.25/m/c. 

Let us note that the photon production rate in (78) can be expressed through 
photon occupation number n{k): 

dR^ 2kT dn{kT coshrj) 

= Wr dTr ^^^^ 

(with two polarizations included). In particular,taking (83) into account one can 
see that if the velocities of the volume elements, as well as proper time interval in 
the first term in rhs of Eq.(78) are small then (78) is reduced to (56) as it should 
be. The estimation of (79) shows that the photon occupation number n{k) in 
Eq.(83) is small numerically, 

n{k) < 1 

That means in particular that the transition radiation is dominated not by the 
photon amplification but by the ground state rearrangement. 

The transition contribution I^l'' in (78) appears essential only at small mo- 
menta kx- So dealing with single-photon distributions one can use (76) for evo- 
lution parameter r{k). The resulting relative strength of the transition radiation 

R^{kT) = 4r^/lQGP (84) 

appears sizable only in the momentum region kx < (20 — 30)Mef independently 
of the time duration of the transition. It is shown at Fig.l. 

Mach better the transition effect is seen in photon correlations (cf (64-65)) 
where it is first order effect with respect to r{k). Let us note that HBT effect 
for photons has now more complicated form than that in (64) because of finite 



time duration |]15| of the process of photon emission from QGP and it will not be 
exposed here. We consider only the transition effect (the third term in (57,64) 
which gives opposite side correlations) estimating its contribution to two-photon 
correlation function in central rapidity region. Suggesting fast transition we can 
evaluate the contribution in the vicinity of fixed proper time r^. So one only has 
to shift to the rest frame of each element of the expanding volume and perform 
77-integration. Then the extention of the invariant correlator in (63) to the case 
of expanding volume takes the form: 

2A;(6±(ki)6±(k2)) = G(kir + k2T)lj:? (85) 

with 

where we neglected n{k) in comparison with unity (see above). Therefore the 
normalized two-photon correlation function is given by (cf (64-65)) 

C(kiT, k2T) = 1 + Chbt + lRl{kT)G\k^T + k^r) (87) 
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with 




8. 



■8 



Iqgp + 




We calculated the ratio i?2(^r) for different transition times St = 0,6t = 
0.5fm/c,6T = 1/m/c up to = 500MeV. The results are shown at Fig. 2. In 
the region kx < lOOMeV the ratio -R2(^t) is sizable for all these 6t being equal 
4.94 at /ct = 0, reaching maximal value i? ~ 6 at /cy ~ 20MeV and falling down 
at kT = lOOMeV to R2 = 1.78 for 5t = 0, /^a = 0.95 for 5t = 0.5/m/c, R2 = 0.55 
for 6t = Ifm/c . At larger transverse momenta the behaviour of the ratio R 
depends strongly on transition time 6t: in interval (200 — 500)Mel^ the ratio 
i?2(^T) rises for 6t = 0, it is approximately constant for 6t = 0.5/m/c and it 
decreases for 6t = 1/m/c. The asymptotic behaviour of the ratio i?2(A;r) at 
large photon momenta kx depends on relationship of the temperature dependent 
single-particle production rate (decreasing like exp(— and the transition time 
dependent evolution parameter (decreasing like exp(— 7rfc5r/2) in our reference 
model). In any case, one can hope to see the effect of the transition in the region 
kx < lOOMev where the peak of the ratio i?2(^T) is always present. 

7 Conclusion 

Estimation of photon emission accompaning transition between quark-gluon and 
hadron states of matter in heavy ion collisions shows that opposite side photon 
correlations can serve as a sign of the transition if transition time is not large. 
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Figure captions 



Fig.l The relative strength of the transition radiation for transition from 
QGP to hadrons. 

Fig. 2 The relative strength of the opposite side two-photon correlations for 
transition time St — 0,dT — 0.5fm/c,ST = Ifm/c (from top to bottom), see 
Eqs.87-88. 
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